#include "PowerFlowCalculation.h"

gsl_matrix* InequalJacobiMatrix(int GenCount,int Node_Num){
    gsl_matrix *res=gsl_matrix_alloc(GenCount+Node_Num*2-1,GenCount+Node_Num);
    for(int i=0;i<GenCount+Node_Num;i++){
        gsl_matrix_set(res,i,i,1);
    }
    
    return res;


}

